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An analysis of the upper critical field in the gauge model is presented. It is shown that a 'strange 
metallic' normal state has implications for the pairing transition even if the resulting superconducting 
state is essentially conventional. The gauge model is considered as one example of an unconventional 
normal state and the optimally doped and overdoped cases are examined. In the optimally doped 
regime B C 2 rises linearly with decreasing temperature to an enhanced value at T = relative to 
BCS. In the overdoped regime the gauge model predicts Fermi-liquid behavior but there is a weakly 
temperature dependent interaction. This can give rise to some upward curvature of the upper critical 
field and an enhanced zero temperature value. 



PACS numbers: 74.20. Mn,74.60.Ec,72.10.-d 

Recent experiments 1 ' 2 on overdoped single layer (es- 
sentially low T c ) cuprate superconductors have made the 
upper critical field accessible throughout its temperature 
range. The results show a number of anomalies when 
compared to the weak coupling BCS theory: B C 2(T) con- 
tinues to rise sharply even at the lowest temperatures 
appearing to intercept the T = axis with a finite slope 
and at a field much larger than would be expected from 
the value of T c . This temperature dependence contrasts 
with BCS theor y whe re B c2 saturates when the magnetic 
length (lg = \/h/eB) is smaller than the electron ther- 
mal length (£/ = hf3vp /2tt). (In this paper a factor of 2tt 
has been included in £/ to simplify formulae.) A second 
order phase transition is always signalled by the diver- 
gence of a susceptibility calculated in the high tempera- 
ture phase at T c , so the temperature dependence of the 
upper critical field is essentially a property of the normal 
state. Thus a non-Fermi liquid normal state might be 
expected to have an unconventional temperature depen- 
dence of the upper critical field. This paper explores the 
consequences on the upper critical field of one such model 
of the normal state — the gauge model. The optimally 
doped region of this model is studied although the sin- 
gular nature of the gauge-field interactions restricts the 
analysis to qualitative conclusions. The overdoped limit 
of this model is also treated to address the above experi- 
ments. The paper is divided into five sections. The semi- 
classical theory as applied to two dimensions is briefly 
reviewed. This defines the terms used throughout the 
paper and enables a discussion of the role of order pa- 
rameter symmetry. The gauge model is introduced in its 
various regimes and the optimal and overdoped regimes 
are treated separately. 



I. CONVENTIONAL THEORY IN 2D 

The theory of the upper critical field has recently been 
reviewed by Rasolt and Tesanovic 3 . The basic theory is 
considered here to address any possible effect of the order 
parameter symmetry. In real space the gap equation near 
T c may be written as 

A(r) = 

\ E jd D r'jd D lV(r - r>)G„ n (I, r)G_ Un (I, r')A(Z). (1) 

An interaction of the form Vkk' = ^2 n v n cosn(8k — 6^) 
for ek and £k' within hu c of the Fermi surfaces is a delta 
function in real space with a width ro « vf/cOc- The 
angular dependence of the interaction acts over distances 
smaller than ro which is smaller than any other length 
scale in the problem. One may therefore write 

A n (r) = A„/ d D lA n (l)K(l,r;P), (2) 

Jl>r 

where K{l,r) = /T 1 £ w „ G w „ (/, r)G_ w „ (/, r) has been 
introduced as the pair propagator. A short distance cut- 
off Tq is also included and the details of the interaction 
on smaller length scales simply renormalize the coupling 
strength to A„ in each angular momentum channel. It 
can therefore be seen explicitly that, neglecting the ef- 
fects of Pauli pair breaking, the symmetry of the order 
parameter has no qualitative effect on the temperature 
dependence of the upper critical field. This is true pro- 
vided Tq is shorter than lg. 

To obtain the well known weak coupling results 4 , one 
makes the semiclassical approximation 5 whereby the ef- 
fect of the magnetic field on the fermion Green's function 
is included by a gauge-dependent Aharanov-Bohm phase. 
Landau level quantization is assumed to be negligible due 
to temperature or lifetime broadening. In the symmet- 
ric gauge this allows the pair propagator in a magnetic 
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held to be written as K (r, r ; B) = exp(ir x r fig )K \\r — 
r'\; 0). In this gauge there is a simple solution 6 to Eq. (2) 
given by A n (r) = A n exp(—r 2 /2l B ) leading to a pairing 
instability when 
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r dr 



(3) 



This demonstrates explicitly that the upper critical field 
is a probe of a normal state correlator. 

In real space the Matsubara Green's function for free 
fermions in 2D is 
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+ ^H^ 2 \r^k' F + 2ik F Lo/v F ) ifw<0, 
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where and are Hankel functions. 
In the limit k F r > 1 the pair propagator is 
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where the Matsubara frequency dependence is included 
as it will be needed later. Putting = one obtains the 
integral equation for the upper critical field 
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dr 



sinh(r/^/) £/ 
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where N(0) is the density of states. Notice that the up- 
per critical field becomes temperature independent for 
Ib <C if- Note also that the equation may be solved 
in both the zero temperature and the zero magnetic 
field limits to yield B C 2(0) = (ir/^)4>o (kBT c /hv F ) 2 where 
7 = exp(0.577) and </>o = h/e, the magnetic flux quan- 
tum. 

In the experiments mentioned above neither of these 
observations appear to hold true. These results depend 
on the Fermi-liquid form of the normal state propagators. 
If the pair condensate is appearing from a non-Fermi liq- 
uid normal state one might expect a new form for the up- 
per critical field. This paper treats the upper critical field 
when there is a pairing transition from the normal state 
of both the overdoped and optimally doped gauge model 
of spin-charge separation. An alternative approach to 
the upper critical field in a spin-charge separated nor- 
mal state has been undertaken by Dias and Wheatley 7 
who assume a form for the normal state Green's function 
based on a ID model with spin-charge separation. 



II. GAUGE MODEL PHYSICS 

In the gauge model 8 the physical electron is decoupled 
into a fermion which carries the electron spin quantum 



number and a boson hole that carries the charge, cJ. CT = 
fr a b r . Underlying strong correlations between electrons 
are reflected in the coupling to a gauge field ('fictitious' in 
the sense that it exists only in this decoupled approach). 
The scalar part of the gauge field, <j), ensures that the 
local density is constant [pf(r)+pt,(r) = 1] and the vector 
potential, a, ensures the local currents cancel [jf(r) + 
jb(f) = 0]. The partition function for the gauge model 
in the presence of an external electromagnetic field A, 
may be written as a path integral over fermion and boson 
coherent states as follows 



Z[A] 



V[a,4>,f,b]e S o 



(7) 



where S[ and Sq are the actions for non-interacting 
fermions and bosons respectively and the interaction is 
given by 

Sint= dr d 2 rj f -a + j b - (a + eA) +i<p(pf + Pb - !)• 



(8) 

This model has three possible 'normal' state phases: a 
strange metal phase ((/^/l^) = 0, (b ) = 0), a spin-gap 
phase ((/jj-f/l^j) 7^ 0, (&o) = 0) and a Fermi-liquid phase 
((/fcf/lfcj.) = 0) (&o) 7^ 0). The superconducting phase 
has both (/^/I^) ^ and (b) ^ 0. Recent work by 
Ubbens and Lee 9 ' 10 suggests that the spin-gap phase is 
thermodynamically stable only for the double-layer com- 
pounds like YBCO. This work does not treat the spin-gap 
phase. Nagaosa and Lee 11 treat the Ginzburg-Landau 
theory of the spin-charge separated system and the their 
findings are compared with those obtained in this paper. 

The electron pair-propagator in the gauge model is 
constructed from the underlying fermions and bosons: 

K e (r,r',p;B) = ^J2 K n( r > r '> # B f) K -n(r', r, (3; B b ). 



(9) 



The electron pair is made up of a fermion pair moving 
from r to r' while a boson pair moves in the opposite 
direction. The boson Matsubara sum ensures that the 
pairs arrive together in imaginary time. The bosons and 
fermions feel different applied magnetic fields due to the 
screening effect of the fictitious gauge field. When calcu- 
lated at a Gaussian level these fields are found to be 12 
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where \b and Xf are the diamagnetic susceptibilities of 
the boson and fermion systems respectively. 
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III. OPTIMALLY DOPED REGIME 



where 



In the optimally doped regime the gauge field is mass- 
less. The convolution of the fermion and boson pair- 
propagators should therefore be performed in the pres- 
ence of the fluctuating gauge field. This is equivalent to 
constructing the electron Green's function in the gauge 
model to derive the pair propagator which has not yet 
been achieved. Instead a form for the upper critical 
field in the optimally doped regime is derived here by 
an approximate consideration of the effect of the gauge- 
field fluctuations. Recent work 13 ' 14 has suggested that 
fermions in a fluctuating gauge field have Fermi-liquid 
response functions at q = 0. It will therefore be assumed 
that the form of the fermion pair-propagator in the pres- 
ence of the gauge field is unchanged from Eq. (5). 

The bosons are strongly effected by the scattering from 
the gauge field. If this were not so, they would condense 
into the lowest Landau level and lead to an upper critical 
field that rises with increasing temperature. As it is, the 
gauge-field fluctuations will wash out the Landau level 
structure making a semi-classical approximation possible 
again. As usual with this problem one can make the 
quasi-static approximation for the gauge- field propagator 
which leads to the problem of bosons moving in a static, 
short-range, random magnetic field. We may obtain the 
boson pair-propagator from the single particle density 
matrix by comparing the density matrix with the Green's 
function in terms of the eigenstates of the problem 



p(r,0) = Y t <l>UO)Mr)e- 0£k 



(ihuj n - € k + H b ) 

p(r,p)e ( ^ +ihu ^dp 



(11) 
(12) 
(13) 



The Laplace transform is well defined for bosons be- 
cause fib < 0. By interchanging the order of summation 
over Matsubara frequencies and performing the Laplace 
transforms, one obtains an expression for the boson pair- 
propagator in terms of the single-particle density matrix 

/>oo °° 

K&(r) = dx^2p(r,x)p(r,x + \m\P) 

J° TO=-00 

x exp[(2/x 6 + iQ)x + \m\fi b P] . (14) 

Since the frequency dependence of the fermion pair- 
propagator is known from Eq. (5), the frequency convo- 
lution of Eq. (9) may be performed yielding the electron 
pair-propagator 
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fits - cos(2irx/P) 



g<yX ^ p Lcosh(r-/^/) + cos(27ra;//3) 



(16) 



In the limit that r <C £/ it can be seen that g(x) ~ 
S(x — n(3). This is the relevant limit when the boson 
length-scale (to be defined later) is smaller than £/ and 
so controls the behavior of the electron pair-propagator. 

To make further progress it is necessary to find a form 
for the single-particle density matrix for bosons in the 
strange metal phase of the gauge model. While this is 
not available, an approximate mapping 15 of the single 
particle problem to the Caldeira-Leggett model of dissi- 
pation 16 may be exploited. This mapping can be sum- 
marized as follows. The problem of particle motion in 
the quasi-static limit of gauge-field fluctuations is char- 
acterized by the Amperean action 
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where So is the free particle action and x = Xf + Xb 
is approximately the Landau diamagnetic susceptibility 
of the fermions (l/127rm/). The best Gaussian action 
that has the same scaling properties on transformations 
of the path may be found by substituting the periodic 
free-particle diffusion law into the log. This leads to an 
action with the same form as that of the Caldeira-Leggett 
model 
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(18) 



Here a general path r T that moves a distance R in imagi- 
nary time has been parameterized as Rt /hp+^2 n r n e %UnT 
where u> n is a boson Matsubara frequency. The form of 
the density matrix for bosons in the gauge model will be 
approximated by the single particle density matrix in the 
above Caldeira-Leggett model. (See Ref. 17 for a critical 
analysis of this approach.) 

The single-particle density matrix in the Caldeira- 
Leggett model in 2D is found by integrating over the 
fluctuations from the straight line path to give 



Pcl(t,P) 



irh 2 p 



x exp 



T(hPuj c /2Tr + l)r (hP/2irT + 1) 

T(hpUJ c /2Tr + hp/27TT + 1) 

~ 2 / 1 ln2Y 



mbr 
~~2T 



HP + 7TT J 



(19) 



Here lo c is a high-frequency cutoff. A simplified den- 
sity matrix which captures the asymptotic forms of the 
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iree-particle limit (hp <C tq) and the dissipative regime 
(h/3 > r ) is 



PCL(r,/3)^-^- + -;]Pxp 



_ ( -Hl±- j- — 

"2 Wp + e b 



(20) 



DifFusion in the dissipative r egime is essen tially cutoff by 
a localization length = -\/7r7iTo/m&ln2. A shift in the 
band edge has been absorbed into the chemical potential. 
The chemical potential is found from the implicit relation 

n b = Y in p(0,n(3)e ni * bP - 

Working in the dissipative regime fifi ^> tq with <C 
£/ the electron pair-propagator is found to be 



N(0)nl exp (-r 2 /$) 
27r£/r sinh(r/£/) 



(21) 



The effective magnetic fields seen by the fermions and 
bosons are included via the semiclassical approximation. 
Since these effective fields of Eq. (10) combine to recon- 
struct the original electron, one may insert the above 
pair propagator directly into Eq. (3) to obtain the upper 
critical field which has the form 



B c2 (T) 



bomi, In 2 
■K 2 h 



[r -\T c )-r -\T)] 



(22) 



Note that the scattering rates that appear are those due 
to scattering from the gauge-field fluctuations and not 
that from impurity scattering which may dominate the 
resistivity at the lowest temperatures. 

The mapping to the Caldeira-Leggett model gives the 
lifetime due to scattering from the gauge field as tq = 
hPmb/3irmf (which leads to a physical resistivity linear 
in T). Thus for reasonable values of the boson to fermion 
mass ratio it is justified to b e working in the dissipative 
regime. Similarly ~ \fk^Tjtp <C 1 so it is clear 

that the boson length-scale controls the electron pair- 
propagator. A form for the upper critical field may be 
now be found 



B c2 {T) = (j>J^k F 



k B (T c -T) 



(23) 



Thus we have a linear increase in upper critical field as 
the temperature is lowered, rising to an enhanced value 
of B C 2(T = 0) compared with that predicted by BCS. 
The zero-temperature upper critical field is enhanced in 
this model by a factor ~ 2kp£o where £o = Tivp/nA, the 
superconducting coherence length. This number is of or- 
der 10 in YBCO for example. There is no saturation of 
the upper critical field within this approach. However, 
this approximation assumes that the bosons never prop- 
erly Bose condense but that ring exchange loops continue 
to grow linearly with fi rather than exponentially. The 
effect of Bose condensation is treated in the next section. 



IV. OVERDOPED REGIME 

The above discussion has necessarily been qualitative 
due to the difficulty of studying bosons interacting with 
a gauge field. These problems are reduced if one works 
at temperatures and doping regimes where the bosons 
have condensed. One might expect this to be the case 
for the overdoped single crystals recently measured. The 
Bose condensate will screen the gauge field so there is 
no long-range interaction present. There are two possi- 
ble mechanisms in this model which could give rise to 
the unusual temperature dependence of the upper crit- 
ical field seen experimentally. One is the effect of the 
convolution of the fermion pair-propagator with the bo- 
son pair-propagator in a strictly two-dimensional Bose 
condensate. The point here is that one normally assumes 
that once the bosons have condensed their field operators 
may be replaced by c-numbers and conventional Fermi- 
liquid behavior is recovered. In 2D one can not make this 
replacement and physical responses may differ from those 
of the Fermi liquid because of the power law correlations 
in the Bose condensate and the frequency convolution 
with this boson response. A second distinct mechanism 
which may lead to an unconventional upper critical field 
is the screening of the gauge-field interaction which ac- 
quires a temperature dependence due to the temperature 
dependence of the effective Bose condensate density. 

The boson pair-propagator in a 2D superfiuid is deter- 
mined using the method of Popov 18 . A boson repulsion 
term u is introduced to precipitate a true 2D superfiuid 
state. Such a term arises naturally in the gauge-model 
description of the t-J model. Popov's method involves 
integrating out the 'fast' fluctuations of the boson field tp 
— the fields which fluctuate above some cutoff wavevec- 
tor ko- This leaves an effective theory expressed in terms 
of the 'slow' fields. The starting action for the bosons is 

r pn r - n 2 - - 

S = / dr d 2 xibd T ib VipVip + pipip 

Jo J 2m 
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(24) 



where the primed fields carry coordinates (x',t) and the 
unprimed fields carry (x,t). The slow fields are conve- 
niently re-expressed in terms of their density fluctuations 
7r [measured from the value at the cutoff p(ko)] and their 
phase variations <p 



1p s l QW (x,T): 
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[tv(x,t) - p(k )]exp[-iip(x,T)] . (25) 



Here b is the boson annihilation operator. The integra- 
tion of the fast variables is performed by Popov and he 
shows that the quadratic part of the effective action of 
the slow variables may be expressed as 
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Here is the boson density and £o is the i-matrix which 
contains the effect of the fast fields to quadratic order on 
the slow fields. One now uses this effective action to de- 
termine the boson pair-propagator in the 2D superfluid. 
One ignores the effect of density fluctuations and retains 
the phase variations to calculate 



K sb (r,r) = (^(r,T)^(r,T)^(0, 0)^(0,0)) 
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The speed of sound for the B ogoliubov spectrum has been 
denoted by v s = ^top s /m b and p s is the superfluid den- 
sity at T = 0. This correlator is the square of Popov's 
result for the Green's function but now generalized to 
include the imaginary time dependence. However one 
can see that in the large distance limit the pair propa- 
gator becomes static. Finally one obtains the result for 
the frequency dependence of the superfluid boson pair- 
propagator 



K$(r)=n 2 6(n) 
= < (r/6) 



irh 2 (3p s 



2h/3v s 

- T/n s{n) 



(29) 



The temperature T b = irh 2 6 jhBm b a^ is the bare boson 
degeneracy temperature which is much larger than T c in 
the overdoped regime. (The number density of bosons 
has been expressed in terms of the doping fraction S 
and the lattice constant do.) A new length scale has en- 
tered £ s = 2h/3v s , the thermal length for long wavelength 
phonons in the 2D superfluid. In order to provide an or- 
der of magnitude estimate for this, note that to is not 
well defined in 2D as k — » so but one may redefine it in 
terms of the chemical potential 19 . An estimate for the ve- 
locity of sound in the dilute limit is v s ss \J Ai:Ti 2 8 / m b a^ 
which, since the fermion and boson masses are of sim- 
ilar order, will be comparable with the Fermi velocity. 
Thus the pair propagator decays with a very weak power 
and is delta functionlike in Matsubara frequency. One 
can conclude that the convolution of the superfluid bo- 
son pair-propagator with the fermions does nothing to 
change the Fermi-liquid form of Eq. (5). 

This shows that including a strictly 2D Bose conden- 
sate, the form of the upper critical field is unchanged from 
the Fermi-liquid result in the gauge model. However this 
omits the residual effect of the screened gauge-field fluc- 
tuations. These fluctuations represent a screened trans- 
verse current-current interaction between the fermions 
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where PT af3 (q) = S a/S — q a q^/q 2 is the transverse projec- 
tor. Unlike the 'strange metal' phase where the current- 
current interaction is singular for low momentum trans- 
fers, here the presence of the Bose condensate gives mass 
to the gauge field and means all momentum transfers are 
important. To understand the effect of this interaction 
on the pairing, consider the scattering of two particles 
on opposite sides of the Fermi surface around an angle 8. 
One finds an interaction 



V kl 



m b e F 
4m fPs (T) 



(1 - cos 261) 



(31) 



which is repulsive in the s-channel and attractive in the 
d-channel and in particular leads to a repulsive pseudo- 
potential p = m b nf /AirifPs for s-wave pairing. By as- 
suming a form for the boson superfluid density one may 
estimate the effect of this temperature dependent inter- 
action on the upper-critical field. A worst case scenario 
is that p s (T) ~ n b (l — T/T b ) which is the intermediate 
temperature form for p s (T) in the dilute limit found by 
Fisher and Hohenberg 19 . Figure (1) shows the conse- 
quences of such a temperature dependent coupling as- 
suming a Fermi velocity of 1.5 x 10 5 m/s, a coupling 
iV(0)A such that T c = 15.35K and a doping of S ~ 0.3. 
An effective boson mass of around 6m/ and T b ~ 5007^ 
gives an enhancement of B C 2(T = 0) appropriate to that 
reported in Ref. 1. This also gives upward curvature of 
B C 2 with decreasing temperature but does not explain 
the dramatic low temperature behavior of B c2 seen by 
Mackenzie et al. 

The above discussion assumes an instantaneous pair- 
ing interaction. If the pairing mechanism is phonon- 
dominated for these overdoped materials then one should 
also worry about the effects of a retarded interaction 
which considerably reduces the effect of the Coulomb 
pseudo-potential in conventional metals 20 ' 21 . As usual 
p* = p /[I + p\og(uj p I lu c )}, where lo p is the frequency cut- 
off scale of the fictitious gauge field. With the bosons 
condensed this frequency scale will be determined by the 
condensate density: lo v « h6/m b aQ. Thus log(w p /a; c ) « 
log(Tf,/l£)) « 2. While the value of the log is small com- 
pared with that for a conventional metal, it is sufficient 
to remove the unconventional temperature dependencies 
of the upper critical field seen in Fig. (1). 



V. CONCLUSION 

In conclusion, the upper critical field has been studied 
here as a probe of the normal state pair-propagator. In- 
dependent of the symmetry of the order parameter, the 
Fermi-liquid result shows saturation of the upper criti- 
cal field at low temperatures. In the gauge model there 
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FIG. 1. The normalized upper critical field in the over- 
doped regime. The BCS result is compared with the effect of 
a temperature-dependent coupling due to scattering from the 
screened gauge field. The Bose condensation temperature is 
30T C and two values of effective mass ratio are shown. The 
value of B C 2(0) is also enhanced relative to BCS by a factor 
of 2 for nib/mf = 1 and 10 for m(,/m/ = 6. 

are at least two possible normal state phases from which 
a superconducting phase may emerge. In the optimally 
and low doped phases where there is no spin gap, the 
upper-critical field is dominated by the short boson ther- 
mal length — a consequence of strong scattering from the 
gauge field. This leads to an enhancement of the zero 
temperature value for the upper critical field — it is found 
to be proportional to T c in contrast to T c 2 for a Fermi 
liquid. The upper critical field in this regime is found to 
have a linear temperature dependence with no saturation 
as long as the bosons remain uncondensed. The slope at 
T c will therefore be enhanced relative to BCS in agree- 
ment with the phenomenological approach of Nagaosa 
and Lee 11 . In the overdoped phase it has been shown that 
the Fermi-liquid temperature dependence should be ex- 
pected (again consistent with the Ginzburg-Landau ap- 
proach). However residual scattering from the screened 
gauge field can lead to some unusual upward curvature 
of the upper critical field and an enhancement of -B C 2(0). 
These effects are almost negligible in the presence of a 
retarded pairing interaction. 
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